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^ . Abstract. By means of two simple examples: phase and amplitude damping, the 

I impact of decoherence on the dynamical Casimir effect is investigated. Even without 

dissipating energy (i.e., pure phase damping), the amount of created particles can be 
0^ , diminished significantly via the coupling to the environment (reservoir theory) inducing 

decoherence. For a simple microscopic model, it is demonstrated that spontaneous 
. decays within the medium generate those problems - Rabi oscillations are far more 

^ ' advantageous in that respect. These findings are particularly relevant in view of a 

recently proposed experimental verification of the dynamical Casimir effect. 

o 

oo 
O 

o 



PACS numbers: 42.50.Lc, 03.65.Yz, 03.70.+k, 42.50.Dv. 
1. Introduction 



pH . Nearly a quarter of a century after the discovery of the static Casimir effect jl], it 
has been reahsed that the (non-inertial) motion of one of the mirrors should induce 
O^' the creation of real particles out of the virtual quantum vacuum fluctuations [2j. 

Unfortunately, this striking prediction (the dynamical Casimir effect) has not yet been 
^ . verified experimentally. 

■ To this end, it should be advantageous to exploit the phenomenon of parametric 

resonance through a periodic perturbation of the discrete eigen-frequencies of a finite 
cavity with the external (perturbation) frequency matching twice the unperturbed eigen- 
frequency of one of the cavity modes, for example. In the exact resonance case, the 
effective (time- averaged) Hamiltonian Hes can be derived by means of the rotating wave 
approximation. If the aspect ratio of the cavity and the resonant mode is chosen such 
that there is no resonant inter-mode coupling (see, e.g., [3j), the effective Hamiltonian 
is just the generator of a squeezing operator {h = c = 1 throughout) 

Hefi = i^[{a^f-a'] , (1) 

with ^ depending on the strength of the perturbation etc., and a\ a denoting the usual 
bosonic creation and annihilation operators for the resonant mode with the well-known 
commutation relation [d, a^] = 1. For example, if the perturbation is induced by an 
oscillation of the dielectric permittivity of a thin dielectric slab with thickness a at one 
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of the walls of the cavity with length L ^ a, the squeezing parameter ^ depends on the 
wave-number fcy parallel to the thin slab and the external (perturbation) frequency lj 
in the following way [Hj 

2 u L 

with X denoting the amplitude of oscillation of the inverse dielectric permittivity. The 
dynamically squeezed vacuum state contains an exponentially growing particle number 

(n) = (0| exp{+iH^st} n expf-ii^cfft} |0) = sinh2(2^t) , (3) 

which should facilitate the measurement (in the ideal case). 

Of course, aiming at an experimental verification of the dynamical Casimir effect 
(cf. [4j), deviations from the ideal behaviour described above have to be taken into 
account as well. Previous studies have been devoted to the effects of detuning (i.e., 
a deviation of the external frequency from the exact resonance) and losses (i.e., a 
finite quality factor of the cavity), see, e.g., 0. However, squeezed states such 
as exp{— iffcffO |0) highly non-classical states (see, e.g., El EI) and therefore 
particularly vulnerable to decoherence in a more general form (see, e.g., [Oj). 

The impact of decoherence on the dynamical Casimir effect will be studied in this 
Article by means of some representative examples for decoherence channels as well as an 
explicit microscopic model. It should be mentioned here that the following investigations 
are devoted to the decoherence of the photon field due to the coupling to the environment 
- not to be confused with the decoherence of the quantised position of a mirror due to 
the photon field (i.e, the dynamical Casimir effect) as discussed in [TUj . 



2. Reservoir 



In the following, we discuss a very simplified model for the reservoir - which nevertheless 
yields some insight into the basic mechanism and illustrates the main idea. Let us 
consider a weak interaction of the mode a\ a under investigation to a reservoir containing 
many degrees of freedom, which will be approximated by independent two-level systems. 
Without any loss of generality, the undisturbed Hamiltonian of the two-level systems 
can be cast (after a suitable internal rotation) into the form X)/ ^i^^i with 2uji denoting 
their intrinsic energy gap and being the usual Pauli spin matrices. The undisturbed 
Hamiltonian of the mode a\ a without squeezing is just VLci^a. In this simple model for 
the reservoir, the most general first-order interaction Hamiltonian reads 

^S=E(«^^/-^/ + «^/-^/) ' (4) 

— * — # 

where A/ denote the small |A/| first-order coupling constants. Since the coupling is 
supposed to be small, we may apply the rotating wave approximation after which only 
the resonant terms (cu/ = fi) survive (see, e.g., [HI EI) 

m E A. {aUl + aal) . (5) 
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Consequently, this effective Hamiltonian consists of hopping operators and describes 
the decay channel. Now let us consider the most general second-order interaction 
Hamiltonian 

= E ((«^)'C/ ■ + a^a A, ■ a, + c?Ci ■ oi) , (6) 

containing the second-order coupling constants C/ and A/ which are usually much smaller 
than the first-order couplings A/. Again applying the rotating wave approximation 

AS E«'«A,a. + E C/ ((a^)Vi + d'ai) , (7) 

we observe that the second term describes resonant two-particle hopping and is usually 
small compared to one-particle hopping H^^^ since the second-order perturbation (j is 
supposed to be smaller than the first order A/. 

However, the first term contains a sum over all reservoir modes - without 
restrictions due to resonance conditions - and thus may well be of comparable magnitude 
or even larger than H^^^. This term corresponds to weak measurements (see, e.g., ^) of 
the particle number d^a without dissipating energy (as [d^d A/cr^, d^d] = 0) and is called 
phase damping. In summary, we obtain the two leading contributions 

= E («^^' + «^+) + E , (8) 

U!j=Q I 

where the first (hopping) term describes amplitude damping (decay channel) and the 
second contribution corresponds to phase damping. 

After averaging over the degrees of freedom of the reservoir by means of the 
usual (e.g., Born-Markov [HE]) approximations, we obtain the master equation in the 
Lindblad form (see, e.g., IS]) 
dg 



—I 
dt 



Aff, ^] + E {^L^fLl - {Lt L„, q]) , (9) 



with two different Lindblad operators La- The first one Li arises from the first term in 
and describes amplitude damping Li = ^/^ d, whereas the second contribution 
in generates phase damping with the Lindblad operator L2 = V^d^d. The two 

damping coefficients 7 and F are determined by the coupling constants A/ and A/ as 
well as the number of involved degrees of freedom of the reservoir. Since this number is 
much larger for phase damping due to the absence of resonance conditions, the phase 
damping rate F can be as large as or even larger than the amplitude damping rate 7 in 
spite of the fact that the second-order coupling constants A/ are usually much smaller 
than the first order A/. 



3. Master Equation and Particle Creation 



For the two decoherence channels discussed above, phase and amplitude damping, we 
shall now investigate the impact of decoherence on the amount of created particles 



Decoherence versus Dynamical Casimir Effect 



4 



h = a'^a. The mode under consideration a^,d is described by the effective density 
matrix g whose evolution is governed by the master equation Q 

j^Q = ^[{a^y - a\ Q\ + - {2nQh - {v? , q\) + | {flaqa} - {n, q\) , (10) 

with the squeezing parameter ^ and the two damping coefficients 7 and F representing 
the amphtude damping rate and the phase damping rate, respectively. 

Instead of deriving the complete solution ^(t) of the above equation, we just focus 
on the most interesting quantity, the time-dependent expectation value of the particle 
number (in the interaction picture with dh/dt = 0) 

±in)^±Trm'^ = Tr[{lm)n]. (11) 

After inserting the master equation (fTUI) . only the first term (^) and the amplitude 
damping term (7) contribute (as phase damping does not dissipate energy), and, after 
some algebra, one obtains 

|(n) = 2^(at2 + a2)-7(n). (12) 

In view of the squeezing generated by Hes, the expectation value (a^^ + a^) will be non- 
zero in general. By using the same method as above, one can derive the time-dependence 
of this quantity 

A(at2 + a') = 8^{n) - (2r + ^){d^^ + d') + 4e (13) 

which fortunately closes the system of equations. This system of two first-order 
differential equations may be solved easily and gives the general solution 

(n)(t) = c+ r + + lee |_^^ ^ ^ yp + ie^} 
+ c^- ^^^^1^ exp {-(7 + r) - ypTI^ 

+ ^ (14) 

217 + 7^ - 16^^ ' 

where C± are integration constants. After incorporating the initial conditions 
(n)(t = 0) = no and ((a^)^ + a^)(t = 0) = (which holds for any initial density matrix 
that is diagonal in the particle number basis such as a thermal ensemble) one arrives at 

4Vr^ + lee r ° v/r^TT6F + r + 7r 



1 16^ 



2 



2r2 + i6^2_ (r + ^)2 • (15) 

Note that the apparent divergence at yJT'^ + 16^ = F + 7 disappears since the 
singularities in the first and last terms cancel each other. 
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As a consistency check, we reproduce Eq. (jH)) for the initial vacuum state no = 
and no decoherence F = 7 = 0. Furthermore, for neghgible phase damping F <^ 
the characteristic exponent is given by 4^ — 7 in accordance with Fe., one only 
obtains exponential growth if the particle creation process (i.e., squeezing) is faster 
than the decay 4^ > 7 (as one might expect, cf. [SI). In the presence of additional phase 
damping, this threshold is shifted to 4^ > 1/7^ + 27F according to Eq. (fT3j) . 

Fet us study the behaviour at very short and late times, respectively. For a positive 
amplitude damping rate 7 > 0, the particle number {n){f) first decreases linearly in time 
(unless no = 0, of course) as the squeezing mechanism sets in at quadratic order only. 
At late times, the particle number either increases exponentially (above the threshold 
4,^ > + 27F) or settles down to a constant (and for ^ > positive) value (below 
the threshold 4^ < ^7^ + 27F). 



4. Decoherence without Dissipation — Phase Damping 



Without amplitude damping 7 = (high quality factor), the situation is qualitatively 
different from the scenario with dissipation. For short times, the number of particles 
increases quadratically and for late times, the particle number is always growing 
exponentially, but the characteristic exponent is reduced. For fast decoherence rates 
F ^ ^ (and a high quality factor), the leading term behaves as 



{n){t) 



no + 



exp 



(16) 



27 ' [ F J 2 

i.e., the particle creation rate is strongly reduced by a factor of order F/^ ^ 1. (Note 
that, for fast decoherence T ^ ^ and intermediate times 1/F -C t <^ T/'C^, one obtains 
quasi-linear behaviour according to the equation above.) 

The reduction of the number of created particles can be understood in the following 
way: The permanent weak measurement of the particle number implies a suppression 
of the off-diagonal elements of the density matrix (i.e., the interference terms) in the 
particle number basis (see, e.g., 0). E.g., without squeezing ^ = 0, we would have 

Qmnit) = {n\ g{t) \m) = Qmn{0) exp |-^("^ - rift^ ■ (17) 

A strong - i.e., decisive - measurement of the particle number would set all off-diagonal 
elements to zero. On the other hand, squeezing is a highly non-classical effect (see, 
e.g., inilZll^) and the related mechanism for particle creation involves the transfer of 
occupation from the diagonal to the off-diagonal elements and back 

^ Qoo Q02 ■■■ \ 



Qr. 



i 

Q2Q 



Qu 



Q22 



V 



(18) 



Therefore, the suppression of the off-diagonal elements diminishes this effect. Roughly 
speaking, we lose the phase information/phase control which is necessary for constructive 
interference and resonance. 
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For fast decoherence F ^ ^, one could visualise the major effect via a gedanken 
experiment combining the time-evolution governed by H^s with permanently repeated 
strong measurements of the particle number at equidistant time intervals At = (9(1/F). 
In the Heisenberg picture, the unitary evolution between two measurements can be 
described by a Bogoliubov transformation 

a{t + At) = a{At)a{t) + (3{At)a\t) , (19) 

with the Bogoliubov coefficients obeying the bosonic unitarity relation |ap — = 1. 
For an initial state g which is diagonal in particle number basis (such as after a strong 
measurement of n), the expectation value evolves as 

{h{t + At)) = {n{t)) + |/5(At)ni + 2{n{t))] . (20) 

In the second term on the right-hand side, one may distinguish the pure quantum 
vacuum effect |/?(At)p and the contribution of classical resonance 2\(3{At)\'^{h(t)). 
For fast decoherence F ^ ^, the time intervals are small ^At ^ 1 facilitating the 
introduction of an approximate differential equation At — »• dt 

|(n(t))cx^[l + 2(n(t))], (21) 

where we have used |/?(At)p = sinh^(2^At) ^ 4^^^^^ g^j^^ ^ C(l/F). We observe 
that the pure quantum vacuum effect d{n{t))/dt oc /T would only yield a linear 
increase oc t^'^/T (permanent reshuffling from diagonal to off-diagonal elements of Qmn 
and accumulating higher diagonal terms), but the classical resonance contributions 
d{n(t))/dt oc 2(n(t))^^/F amplify all particles present and thereby lead to an exponential 
growth with a reduced exponent oc t^^/F as in Eq. (fTI)|) . 



5. Microscopic Model 

In order to apply the above results to an explicit scenario, one has to estimate the 
magnitude of the involved quantities - such as the phase and amplitude damping rates 
F and 7, respectively. This can be achieved either by experimental means, i.e., measuring 
the quantum coherence time (oc 1/F) and the quality factor (oc I/7) of the used cavity 
etc., or via theoretical calculations - based on a microscopic model, for example. In the 
following, we shall present a very simple microscopic model for a material with a time- 
dependent index of refraction - which, nevertheless, yields some interesting conclusions. 

Let us consider a specifically designed semi-conductor at low temperatures, for 
example, whose localised single valence electrons (far below the conducting band) 
can occupy three levels described by the amplitudes 4'a,b,c(t)- AH the other levels 
are supposed to be unimportant and are not occupied at all. Illuminating the semi- 
conductor with a strong external Laser beam tuned to the frequency of transition from 
the lowest-lying electronic state a to the first excited state b, we can manipulate the 
amplitudes of the electronic states ipa,b{t). In addition, a small and slowly varying 
electric test field £ (e.g., micro- waves) is acting on the three- level system. The spatial 
distribution of the ground state a is supposed to be very compact such that it is basically 
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unaffected by the electric test field £{t)\ but the excited states h and c are more spread 
out and hence couple to £{t) with the dipole moment k (dipole approximation). In the 
rotating wave approximation, the Lagrangian of the described three-level system reads 

C = irai^a + iVbi^b + - Ao; rci^, + (/« ^(t)^,*^c + ^KitWai^b + h.C.) , (22) 

with Ao; denoting the energy gap between states h and c, and VL^{t) the Rabi frequency 
of the strong external Laser field (see, e.g., [3 IE])- Since £{t) is small, we may use linear 
response theory {ip^ <^ 4'a,i^b), and employing the adiabatic approximation* 

KS{t)'il]b = Auj ipc — ii^c ^ ''Pc , (23) 

because £{t) is supposed to be slowly varying, we finally (after ehminating ipc) arrive at 
the contribution to the effective Lagrangian for the test field £ 

^eff = mt)?^S\t) = '-^^^^£\t). (24) 
Au; 2 

Consequently, by means of a Laser-induced excitation of the polarisable 6-level ipb{t), 
we may generate a time- dependent effective dielectric permittivity Ees- This, in turn, 
can serve as a perturbation for the dynamical Casimir effect |31 E] • 

For a rough estimate of the maximum possible order of magnitude of SeS, we assume 
one three-level system per lattice site with a dipole couphng k corresponding to the 
lattice spacing of a few Angstroms (i.e., the electronic states b and c are spread out over 
a few Angstroms). The Rabi frequency f^R determines the time-scale for the dynamics 
of the micro-wave test field £ and should be around a few GHz in order to achieve 
resonance The adiabatic approximation employed above requires the energy gap 
Auj to be much larger than that frequency. If we tune the three-level system in order to 
make this energy gap as small as possible within the region of validity of the adiabatic 
approximation, say (9(100 GHz), the effective dielectric permittivity ^cfr can be as large 
as Ees < O(IO^), which is more than sufficient. 

6. Spontaneous decay noise/decoherence 

So far, we only considered stimulated transitions a ^ b induced by the (external) Laser 
field. Under this assumption, one obtains an effective Lagrangian ()24|) leading to a 
unitary evolution, i.e., no decoherence. As it will become evident below, this result 
changes drastically for a spontaneous decay back to the ground state a instead of a 
stimulated transition. 

Assuming that the spontaneous decay (involving optical frequencies) occurs much 
faster than the (slow) dynamics (e.g., micro- wave frequencies) of the test field £{t), we 
may adopt the sudden approximation and omit the time-dependence of £(t) during that 
process. The influence of the field £ perturbs the polarisable 6-state inducing a mixing 

* Outside the region of validity of the adiabatic approximation, Alo should be replaced by the detuning 
between the frequency of the test field £ and the b c transition frequency Aw. 
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with the c-state and thus an energy-shift, which can be calculated using second-order 
stationary perturbation theory 

=Er + (61^16) + (25) 
-fcfe - Ec 

The first-order term vanishes in the dipole approximation (cf. the selection rules) 
{h\£\h) = 0. Depending on the slowly varying electric field S{t), the spontaneous 
decay of the excited state - which is a mixture of the states b and c - releases a different 
energy El [£{t)] than the originally absorbed Laser photon in general. Note that the 
totahsed energy shift E^^^ — Ef^ in Eq. exactly corresponds to the polarisation 
term in Eq. ()24p . i.e., the polarisation energy £V of the medium 

eeff^^ =£V = £'^ + £V. (26) 

So far, the electric field £ has been treated as a classical perturbation in Eq. ()25j] which 
is strictly valid in the classical limit (coherent state with many micro-wave photons) 
only. For the quantised micro-wave radiation field in a dynamically squeezed state, the 
situation is more complicated. Expressed in terms of the relevant mode a\ a with the 
frequency f2, the electric field operator can be estimated by 

£^ = o(^\a^ + af] , (27) 



y 

with V denoting the volume of the cavity. In a squeezed state, the expectation 
value of £'^ and its variance are strongly time-dependent in view of a){t) = e*'^*a'l"(0). 
Since we may decompose the squeezed state into a linear superposition of (nearly 
classical) coherent states, the energy of the emitted photons will fiuctuate randomly 
(depending on the moment of decay). Therefore, spontaneous decays inevitably lead to 
decoherence - first, by generating energy noise (total energy conservation), and, second, 
by effectively performing a weak measurement, since the emitted photon carries away 
some information about the electric field (encoded in its energy, for example - see the 
discussion below). Comparing the energy shift in Eqs. ()25|1 and with the squeezing 
parameter in Eq. we may estimate the total noise energy per cycle 

E^oise = OiNO , (28) 

with being the total number of micro-wave photons inside the cavity. This noise 
energy is of the same order of magnitude (with the explicit pre-factor depending on the 
details of the geometry and the dynamics of Ees etc.) as the energy gain due to the 
dynamical Casimir effect in each cycle - which already demonstrates one of the main 
drawbacks of spontaneous decays. 

The estimate of the phase damping rate T generated by spontaneous decays is 
rather complicated. To this end, it is necessary to study the entanglement between the 
micro-wave mode a\ a under consideration and the reservoir, cf. p. One indicator is the 
amount of information about the quantum state of the micro-wave field (in particular 
the particle number) leaking out - which can be used as a lower bound for Fdecay Ergo, 
the question is: how much can we learn about the number of micro-wave photons inside 
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the cavity by looking at the (optical) photons emitted during the spontaneous decays? 
Since the spontaneous decays have to be repeated with a frequency of order VL in order to 
achieve resonance 0, the line-width will also be of order VL. Let M denote the number 
of absorbed Laser photons - and hence excited three-level systems. Of course, the 
same number M of (optical) photons must be emitted by (approximately independent) 
spontaneous decays in each cycle and hence the variance of the total emitted energy is 
0{\fMVL). This quantity can be contrasted with the accuracy necessary for measuring 
the total number of micro-wave photons inside the cavity (by looking at the optical 
photons emitted). In the presence of micro- wave photons with an energy of order 
Vt (resonance condition), the total energy difference in Eqs. (f^H|) and is of order 
A^.^, cf. Eqs. (j2I) and ()27|) . The ratio of these quantities N^/{\/Mil) determines the 
information gained in each cycle (i.e., spontaneous decay). As a result, this special type 
of decoherence becomes important if the number of micro-wave photons A^ exceeds the 
threshold set by the number M of absorbed Laser photons 



However, as demonstrated in Section|2l there are far more possibilities for inducing phase 
damping. For example, if the electrons are pumped (by the Laser) into the conducting 
band (in contrast to the three- level system discussed above), many more degrees of 
freedom are available for carrying away energy and information. 

The undesirable consequences of spontaneous decays can be avoided by inducing 
stimulated transitions which are much faster than the spontaneous decay rate - for 
stimulated transitions, the emitted photons have the same quantum numbers as the 
incident light and thus do not carry away energy or information. E.g., for a Laser beam 
with a constant intensity Qb. = const, we obtain the well-known Rabi oscillation 



Rabi oscillations (i.e., stimulated controlled transitions a ^ b instead of a spontaneous 
decay b,c a) maintain coherence - e.g., they are used to manipulate QuBits IHj. 

7. Conclusions 

The objective was to study the impact of decoherence on the dynamical Casimir effect. 
For the pure decay channel (amplitude damping), the dissipation of energy results in a 
subtractive reduction of the characteristic exponent, i.e., the number of particles only 
grows exponentially if their creation (via squeezing) is faster than their decay (given 
by the quality factor of the cavity), cf. However, even without energy loss, e.g., 
for pure phase damping, the dynamical Casimir effect is diminished by decoherence - 
but in a different way. For pure phase damping, the number of created particles always 
grows exponentially, but with a reduced exponent. Furthermore, the pure quantum 
vacuum contribution (i.e., the original dynamical Casimir effect) would only yield a 
linear increase of (h) ; the (reduced) exponential growth is caused by classical resonance 
amplifying already present particles. For an experimental verification of the dynamical 




(29) 



(30) 
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Casimir effect, a large enough quality factor of the cavity is necessary, but not sufficient^ 
the quantum coherence time is also a very important quantity. 

In order to apply the above results to a concrete example, we considered a three- 
level system as a simple microscopic model allowing the coherent control of eeff(^) via 
Laser illumination (of a semi-conductor, for example). As long as the Laser beam can be 
treated as a classical external field (stimulated emission only), the evolution is unitary, 
i.e., without decoherence (under the assumptions made). Spontaneous decays, on the 
other hand, inevitably generate decoherence: firstly, by dissipating the polarisation 
energy of the medium in Eq. (j^Bj) : and, secondly, by allowing information to leak out - 
which effectively corresponds to weak measurements. Although the explicit decoherence 
rate depends on the concrete realisation, spontaneous decays are therefore not desirable 
for a controlled experimental verification of the dynamical Casimir effect. 

As one possible solution for this problem, one could use Rabi oscillations which 
are faster than the rate of spontaneous decays and automatically generate a harmonic 
dependence of Sesif), cf. Eq. (jHUjl . maintaining coherence. Another advantage of Rabi 
oscillations lies in the fact that they do not deposit energy into the material, i.e., there 
is (ideally) no heating - avoiding unwanted excitations which generate decoherence. 

If one departs from the three-level system consisting of localised electrons and 
excites electrons in the conducting band, similar difficulties arise. Apart from 
the aforementioned problems owing to spontaneous decays, a finite conductivity 
always entails additional noise and decoherence according to the fluctuation-dissipation 
theorem. Furthermore, one should be very careful not to excite degrees of freedom 
generating a non-linearity which induces a coupling to higher harmonics. E.g., via 
parametric down-conversion, such a coupling might overcome the frequency separation 
between the Laser and the micro-wave photons - and thereby completely swamp the 
desired signal. 
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